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Abstract: Distance has been one of the basic factors in manufacturing and control fields, and 
ultrasonic distance sensors have been widely used as a low-cost measuring tool. However, 
the propagation of ultrasonic waves is greatly affected by environmental factors such as 
temperature, humidity and atmospheric pressure. In order to solve the problem of inaccurate 
measurement, which is significant within industry, this paper presents a novel ultrasonic 
distance sensor model using networked error correction (NEC) trained on experimental 
data. This is more accurate than other existing approaches because it uses information from 
indirect association with neighboring sensors, which has not been considered before. The 
NEC technique, focusing on optimization of the relationship of the topological structure 
of sensor arrays, is implemented for the compensation of erroneous measurements caused 
by the environment. We apply the maximum likelihood method to determine the optimal 
fusion data set and use a neighbor discovery algorithm to identify neighbor nodes at the top 
speed. Furthermore, we adopt the NEC optimization algorithm, which takes full advantage 
of the correlation coefficients for neighbor sensors. The experimental results demonstrate 
that the ranging errors of the NEC system are within 2.20%; furthermore, the mean absolute 
percentage error is reduced to 0.01% after three iterations of this method, which means 
that the proposed method performs extremely well. The optimized method of distance 
measurement we propose, with the capability of NEC, would bring a significant advantage 
for intelligent industrial automation. 
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1. Introduction 

Distance is one of the most basic factors in manufacturing and control fields. It is used for local 
positioning, object identification, automation control, human-computer interaction, and so on [1,2]. 
In order to measure distance, ultrasonic distance sensors (UDS) have been widely used as a low-cost 
solution [3], and large numbers of them have been set up intensively as sensor arrays [4]. For example, 
in the production lines of assembly workshops — to ensure the size of products, and in parking guidance 
systems for shopping malls — to detect and classify automobiles. However, the propagation of ultrasonic 
waves is greatly affected by the environmental factors such as temperature, humidity, and atmospheric 
pressure [5]. Many traditional methods [6] have been employed to compensate for errors introduced by 
the sensor design. Furthermore, with the development of intelligent industrial automation, sensors have 
been networked to interchange data, and can easily gain supplementary information from neighboring 
nodes. In combination with these features, this paper proposed a novel design of sensor arrays that can 
correct its own errors with information from the network. 

Figure 1. Scenario of production lines, (a) Parallel lines in a workshop, (b) One line in 
a workshop. 




(a) (b) 



As shown in Figure 1 , taking the scenario of production lines as an example, there are many parallel 
lines in a workshop. A number of UDS are set up in each line to measure the size of the semi-finished 
products, and the measured values are passed to a data center (DC) to judge the quality through 
the network. The measured values are unavoidably imprecise, because these sensors are in different 
micro-environments. To solve this issue, each UDS can be equipped with the temperature and humidity 
sensors to correct errors caused by the environment. However, this requires a greater hardware cost, 
more redundant circuit designs, and more complex compensation algorithms. Nevertheless, with the 
information-based development of modern industry, the data from the sensors can be gathered by the DC 
for centralized interactive control. A hardware upgrade of UDS can be avoided — decreasing the total 
cost — if an effective method of DC data analysis and subsequent error correction is carried out. 
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The most basic method of UDS is using a simple pulse echo sensor, which acts as both transmitter 
and receiver. An echo is detected by the receiver after the transmitted pulse propagates outwards. The 
distance between the transmitter and receiver is obtained by counting the elapsed time from the start of 
the transmission to the end of the receipt, defined as the time-of-flight (TOF) [7]. Given that ultrasonic 
speed is a fixed value in a certain environment, the distance can be calculated by the following Equation: 



where d s is the measured distance, v s is the ultrasonic propagation velocity in the medium, and AT is 
the signal TOF. Then, the precision of d s depends on the precision of v s and AT. In particular, the 
distance between the sensor and the floor is a fixed value d 0 if nothing is measured. Therefore, there is a 
zero mode (Z-MODE) in the device meaning that, irrespective of the measured value from the UDS, the 
distance is regarded as d 0 without any definite errors. 

The fact that the industrial environments of sensors are non-ideal implies that it is sometimes difficult 
to detect the real size of the semi-manufactured products. In order to improve the detection process, 
many mathematical methods have been used, closely followed by the use of UDS growing in an 
exponential manner since the 1980s [8]. On one hand, using a comparative and collaborative approach 
is beneficial to multi-analysis and produces compelling results. In [9], it not only presents an optimized 
ultrasound-driven method of incorporation of the amplitude modulation and phase modulation of the 
transmit-receive technique, but also proposes two optional algorithms to achieve TOF calculations, 
according to different signal-to-noise ratio situations — named the zero-crossing tracking and time-shifted 
superposition methods. In [10], it applies the discrete extended Kalman filter to the acquired ultrasonic 
signal, in order to accurately estimate the shape factors of the echo envelope as well as to locate its 
onset; it is also possible to ensure reduced bias and uncertainty in critical TOF measurements. In [11], 
it introduces the received ultrasonic wave peak time sequences of two slightly different frequencies, and 
the experimental results demonstrate that the relative TOF can be accurately identified with much better 
resolution than the wave period. On the other hand, modeling based on the network topology of sensor 
arrays has already been applied in some modern industrial cases. 

As mentioned above, in some of modern industrial cases , these devices have formed a network with 
wired or wireless methods passing data to a DC. Furthermore, many data mining and mathematical 
and statistical algorithms have been used for data analysis. A model using data mining algorithms 
which then processes data based on the same model has been widely employed in industrial fields to 
gain accurate measurements. In [12], it uses the statistics of the relation of selection and theoretical 
variances, the simultaneous operative testing of mathematical expectation, and the variance of errors 
in the one-dimensional measurements. Then in [13], an approach to characterize the model prediction 
errors using a Gaussian mixture model is proposed, which mainly consists of choosing the relevant 
variables to form the error data, optimizing the number of Gaussian components required for the error 
data modeling, and fitting the Gaussian mixture parameters with an expectation-maximization algorithm. 
In [14], it proposes a two-mode adaptive fuzzy controller with an approximation error estimator. In 
the learning mode, the controller employs some modified adaptive laws to tune the fuzzy system 
parameters and an approximation error estimator to compensate for the inherent approximation error; in 
the operating mode, the fuzzy system parameters are fixed — only the estimator is updated online. These 
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methods have been effectively used for data analysis, so they can also be incorporated into networked 
UDS systems for error correction. In general terms, previous studies have suggested that using methods 
that focus on optimization of a single UDS device can improve the performance of distance measurement 
to some extent, however, contradictions among the accuracy, hardware complexity, data update rate and 
measurable range still exist. 

In this study, based on the fact that UDS devices have been organized as a network to gather real-time 
data into a DC, we focus on the mathematical methods to calculate the precise correction value by 
analysis of the data correlation coefficients via networked error correction (NEC). The main features of 
the method are briefly described in this paragraph, where the model of NEC and the analysis of the data 
correlation coefficients are also presented. Field tests are implemented in this study, because they are 
more reliable and more likely to reveal potential problems. To the best of our knowledge, we identify 
for the first time that a method using network topology information and leveraging recent advances in 
optimization would be a more useful technology for distance measurement. 

Over all, the contributions of this paper are summarized as the followings: 

• We have modeled the NEC system on the optimization of the relationship of the topological 
structure of the sensor arrays, and developed a Two-Step Error Correction Process (TSECP) based 
on the relations between neighboring nodes. 

- An abstract model is set up to illustrate the ability of the compensation for error 
measurements caused by the environments. 

- TSECP is described with the theory of fusion set. To compare with the previous NEC 
methods, we have analyzed the design comparisons and the system improvements. 

• We have designed the UDS with the maximum likelihood method, neighbor discovery discovery 
algorithm and the NEC optimization algorithm. 

- We apply the maximum likelihood method to determine the optimal fusion data set and use 
a neighbor discovery algorithm to identify neighbor nodes at the top speed. 

- We adopt the NEC optimization algorithm, take full advantage of the correlation coefficients 
for neighbor sensors. 

The experimental results demonstrate that the NEC system has excellent anti-interference 
performance, and its ranging errors are within 2.2%. Furthermore, the mean absolute percentage error is 
reduced to 0.01% after three iterations of this method. 

Because we take advantage of data correlation coefficients for the first time, the novel NEC model is 
more accurate than existing approaches. Its superiority over the former methods is obvious: instead of 
the solo node, all of the topological information of the sensor arrays is involved in error compensation. 
Observational evidence has frequently been linked to numerical error and propagation of the ultrasonic 
wave affected by the environment; however, there is little direct trial evidence to train and verify the 
exact equation of the two parameters. In view of the significant potential value of this issue, further 
investigation of the equation is warranted. 
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2. Related Work and Motivation 

2.1. Distance Measurement by Ultrasonic Sensor Arrays 

Recently, distance measurement based on UDS has increased significantly. In [15], according to 
the multi-sensor information fusion theory, it uses Kalman filter principle to deal with single ultrasonic 
sensor information and multi-ultrasonic sensor information fusion. The results show that multi-sensor 
distance measurement is more accurate and stable than single sensor distance measurement. In [16], 
it also proposes the experimental dimensionless curves/corresponding equations. Then, it presents the 
relative profile of calibrating coefficients and ambient temperature of the specific ultrasonic sensor to 
accelerate the field calibrating process. Obviously, it is unavailable before the formal stage measuring, 
so the calibrating accuracy with equation estimation is slightly inferior to this process. 

Moveover, the performance of the UDS arrays system can be greatly improved by adopting data 
fusion technologies and self-configuring scheduling protocols. In [17], it presents a self-configuring 
protocol for time- slot allocation problem, under the constraint that no explicit communication is allowed 
in the sensor system. Simulations with two sensor geometries have demonstrated the effectiveness of 
the proposed protocol. In [18], it presents a novel approach toward the identification of odors/gases 
using game theory for feature selection. Considering that the merit of this approach lies in the fact that 
each dimension of the raw data gets adequate representation in the form of Shapley value assignment 
based upon its contribution to class separability. However, it can not be compared to a standard sample 
considering the lack of appropriate scaling of test samples. 

Notably, the implementations of UDS arrays have obtained increasingly influence in the industrial 
world. In [19], based on the two-frequency continuous-wave phase-shift method, it has developed a 
high-resolution UDS system using vernier caliper phase meter, which can reduce the cost and complexity 
of the circuit elements, as well as decrease the noise disturbance induced from the high-frequency clock 
signal. In [20], it presents a localization system for Wireless Sensor Networks (WSN) based on ultrasonic 
TOF measurements. The mean error and standard deviation of the angle and distance measurements lead 
to a small mean localization error of 4.21 cm and a standard deviation of 0.57 cm. But there are also 
some disadvantages: ultrasonic consumes more energy, it has a lower range; and the speed of sound 
depends on the environmental temperature and humidity. 

2.2. Networked Error Correction with Mathematical Statistics 

NEC is the key stone for the implementation of UDS. In [21], it introduces the concept of 
co-integration, describing an equilibrium relationship among a set of time-varying variables, which can 
be utilized to represent the co-integrated relationship through an error correction model and examine 
the effect or local influence on the error correction variable. From another perspective, methods based 
on the topological structure relationship of components have been applied to data analysis. In [22], a 
forward error correction scheme for WSN has been developed to avoid retransmission. It not only saves 
energy but also extends the functionality and enables to handle "BURST ERRORS". It is more effective 
because the energy expended for the retransmission in the case of the original correction scheme is 
saved when the new proposed error correction scheme is implemented. In [23], in order to address the 
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fact that distance measurements are often corrupted by large errors which will lead to inaccurate position 
estimates without correction, it has proposed a scheme of joint node localization and error correction for 
the distance measurements of WSN. Compared with the centralized algorithms, decentralized methods 
are more robust and scalable. However, it does not studies the case that the position information of 
anchors are inaccurate. In [24], it has proposed an adaptive forward error correction algorithm for the 
best effort WSN. This solution does not require additional feedback channels, which makes it compatible 
with WSN. However, it has not justified the choice of link quality metric for this work and gained a 
deeper understanding of the contribution. 

Moreover, mathematical statistics has been widely used for NEC. In [9], it presents a two-level 
neural network fed with the acquired samples of each echo received, employed for post-processing. 
A competitive layer has been used to determine which period of the echo waveform made the threshold 
detector fire; then, a multi-layer perception used to determine the actual echo position. These methods 
focus on the optimizations in the single UDS devices, which put extra pressures on the computation of 
the source-limited UDS device. In [25], it has proposed two estimators based on correlation parameters 
for the two key steps of a practical distributed source coding scheme: first, the computation of the 
side-information at the receiver side; second, the estimation of the required number of bits to compress 
the readings in order to guarantee a certain symbol error probability. This method allows to reduce the 
number of transmitted bits and hence reduces the energy consumption. In [26], it firstly studies the error 
between the value obtained from the sample data of all or a part of nodes, and carries out algorithms 
to assign confidence levels to child nodes of the root on the basis of their corresponding confidence 
levels. The experimental results illustrate that this scheme can save energy extremely. However, it has 
not analyzed the error and confidence level allocation when the beginning time of each node to sample 
data is asynchronous. 

2.3. Motivation 

The observations above demonstrate that measured values are unavoidably imprecise, because of 
the instability of the environment. This motivates us to propose an approach of minimizing the error 
between the measured response and the desired response. This study provides a potential solution based 
on NEC for the reliable and low-cost distance-measurement applications by using the information from 
neighboring sensors — which has not been considered before. This also has the potential to shed light on 
the automation fields. 

3. Process of Networked Error Correction 

3.1. Basic Architecture 

The basic architecture of networked sensors is shown in Figure 2. Each node has a fixed position 
and measures the distance independently. All of them constantly pass the measured data to the DC. 
Besides positional information, two other aspects are considered with changes in time. First, at some 
time point t all of the nodes measure the raw data without any correction hardware (temperature or 
humidity sensors). If a node is in Z-MODE, its value can be compared to the precise value do, and the 
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error of the node can be found by a simple subtraction. Second, in a time interval, t G [ti,^]* eacn 
node measures an independent series of data points (because of the different objects measured). There 
is no distinctive function for a node; however, for the nodes in their fixed positions, the errors caused 
by the environment have the implicit function. In mathematics, an implicit equation is a relation of the 
form R(x\, x n ) = 0, where R is a function of several variables ( often a polynomial), for example, 
x 2 + y 2 — 1 = 0. While, a distinctive function is a relation of the form y = f(x), where y can be 
described directly by equation containing independent variables, for example, y = sin(x). Therefore, if 
a value is corrected, other values may be corrected accordingly. 

Figure 2. Block diagram of the basic architecture. 
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For example, if Nodes / and J are two separate nodes. At some time point, t, they measure the 
raw values I t and J t . If no object has passed through Node /, I t can be compared to d 0 — the error 
in this point is E It = (d 0 — I t ). Further, if over some time interval, Node / and Node J measure at 
three times points t\, t 2 and t 3 obtaining two sets of raw data, {Ii, I 2 , 1 3 } and { J\, J 2 , J3} with errors of 
{£71, E I2 , Ei 3 } and {Eji, E J2 , E J3 }, (Err(I) and Err(J)). Although Ii has no direct relation with I 2 
and J 1, the environments are very similar for both of nodes at the time points, t\, t 2 and t 3 . The errors 
caused by the environment have the inner relations that can be concluded from the trends of Err (I) and 
Err(J) over time, i.e., Err(J) = Fun(Err(I)). If I\ is in Z-MODE, En can be calculated directly, 
then Eji can be calculated from the function, and J x can then be found. 



3.2. Models of NEC 

In the air with a constant pressure, v s (shown in Equation 1) is given by the equation [27]: 



^ (2) 
M v ' 

where R is the gas universal constant, 7 is the specific heat ratio — the ratio of gas pressure fixing specific 

heat capacity in specific heat at constant volume, and M is the molar mass of the gas. The three variables 

all have minute changes in common [28]. T is the thermodynamic temperature of the gas (in Kelvin), 

which can be calculated from the temperature in centigrade t c by T = cq + t c . cq is the Absolute 

Temperature, a constant. Using Equation 2, it can be seen that the ultrasonic velocity will increase by 

approximately 0.6 m/s if the temperature increases every 1 °C. If the default velocity is around 340 m/s 

at 20 °C it is then the primary origin of measurement error in the UDS. 
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In addition, v s is also influenced by other factors of the environment, such as humidity and 
atmospheric pressure. As for AT, it also has unavoidable errors, arising from the attenuation of the 
ultrasonic waves and the features of the detected objects. 

If v s has the error of o(v s ) and AT has the error of o(AT), then the corrected d s (as d s ) can be 
given by: 

i _ (vs + ojvs)) x (AT + o(AT)) 

(3) 

_ v s xAT o(AT) o(y a ) o(y 8 ) x o(AT) 

2 2 s 2 2 

where, o(AT) x v s and o(v s ) x AT are relations; it can be represented as: 

d s = d s + err(a, (3, 5) (4) 

in which d s = and err (a, /3, 5) = ^ xt) s + ^xAT+ °M*°(Ar) = a + (3 + 5. Considering 

that both o(v s ) and o(AT) are infinitesimally small, and the infinitesimal order of 6 is lower than a 
and /3, the value of err(a, (3, 5) heavily relies on a and (3. Equation 3 also demonstrates that there is 
heavy coupling between v s and AT in both a and (3. In order to decouple v s and AT, we decided to train 
the equation of err (a, (3, 5) through sufficient experimental data fitting, using statistical algorithms. 

Mathematical statistics has been widely used in the collection, analysis, interpretation and 
presentation of data in the process of error correction. Variance is used as a measure of discreteness. 
A set of numbers has a probability distribution, expressing how far the numbers deviate from E [X\. 
Variance is a parameter to describe either the actual probability distribution of an observed population 
of numbers or the theoretical probability distribution of a sample (a not-fully-observed population) of 
numbers. If a variable X has the expected value (mean), p = E [X], then the variance of X is given by: 

a 2 = Var(X) = E [(X - [if] = E [X 2 ] - (E [X]) 2 (5) 
where a 2 is the variance. 

For two variables, covariance is a measure of the probability that they change together. The sign of 
the covariance shows the tendency in the linear relationship between the variables, and the magnitude of 
the covariance shows the strength of the linear relation. We define the covariance of the two variances 
X and Fas Cov(X,Y): 

Cov(X, Y) = E[(X - E[X])(Y - E[Y})} = E[XY] - E[X]E[Y] (6) 

While the magnitude of Cov(X, Y) is dependent on the magnitude of X and Y, it cannot compare 
the linear relationships between different variables, so the Pearson correlation coefficient (p) is defined 
to measure the correlation (linear dependence) between two variables, X and Y, giving a value on the 
interval [—1, 1]. If Var(X) ^ 0 and Var(Y) ^ 0, p is defined as: 

p{X ,Y) = , °f X >p (7) 
^Var{X)^Var(Y) 

Considering that p(X, Y) represents the relationship between X and Y, p can be used to quantize the 
relationship between an arbitrary node and its neighboring nodes. Data of some nodes can be corrected 
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by data of others after p(u, v) has been obtained, because the relationship of two nodes u and v, chosen 
at random, is known. 

An abstract model is set up to illustrate the problem. Assume all of the UDS distributed in a large 
area, form a matrix and are at time t; each of them can measure a value, a^, in the position of i and j; 
then, A is the raw value matrix 



Old; 



an, 



a ln 



(8) 



Some of these values could have been corrected by other methods, such as the Z-MODE of a UDS 
device. Assume that aj~j is the corrected value in the position of % and j, then B' is the corrected value 
array. For example, if a 00 , a 01 , an, a mn are corrected to a^, <%j~, a^, then B' has four elements, 
B' = {&o, &i, 6 2 , ^3} = {ooo, ooi, oil, It can also be expressed as a sparse matrix relation: in a 

positron matrix 5 and value matrix is B: 



S t = 



B t = 



' 1, 1, 

0, 1, 

. 0, 0, 

a oo> a 01> 

aio, an, 



0 
0 



(9) 



°0n 
O In 



(10) 



If P = 1 (each value in the matrix is one), then B = B 0 ; it is an initial state. Using simple subtraction, 
the error matrix is easily obtained: 



Err t 



a 00 ~ a 00> a 01 — a 01) 

0, an — an, 



0. 



0. 



0 
0 



0 5 ®"mn Oir, 



(ID 



In this paper, the key information is that a^ is correlated with its neighboring nodes, which may be 
'one-to-one" or "one-to-many" relations. We define the neighbor array: 



(12) 



N = {no,rai,ra 2 , • } 

= {fl(i-i)(,-l), °(0(j-i)' °W(j+i)' ' ' ' ' a (i+l)(i+l)} 

For ay, if a correct function, Corij(N), has been found from the analysis of the history data, then the 
corrected (a^) is 



dij = Cor i:j (N) 



(13) 
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Note that each aij has its own correction function, Cor^N), and the function may also be 
dynamically corrected along with the accumulation of the history data. After the correction, the matrix 
B will be converted to a corrected matrix C: 



a 0(b a 01) 

aio, an, 



0>ln 



(14) 



According to the analysis above, the key problem is to find the suitable Cor^N) for each a^. 
3.3. Relations between Neighboring Nodes 

Neighboring nodes are located in fixed static positions and have similar trends with environment 
changes, so they may have some relations. As shown in Equation 7, p is defined as a measure of the 
correlation, which gives an indication of the strength of the linear relationship between the two random 
variables. If p = 0, then X and Y are uncorrected, and the stronger the correlation, the closer \p\ is 
to one. 

As a simple case, in a time interval, [ti, ^2], where each node of n nodes has collected 1,000 raw 
measured values. As mentioned above, the mathematical expectation matrix of the rough data of each 
node is Ei 



E; 



a 0 



(15) 



where i = 0, 1, 2, .... 



From Equations 6 and 7, the Pearson Correlation Coefficient (p) matrix is defined as p 



Pi,j 



1 


Po,i 


P0,2 • 


' ' P0,n-1 


Pl,0 


1 


Pl,2 


' ' Pl,n-1 


P2,0 


P2,l 


1 


' ' P2,n-1 



(16) 



Pn-1,0 Pn-1,1 Pn-1,2 
where i, j = 0, 1, 2, ... 

Of particular note is that the sum of each line in the p i: j matrix is not likely to be 1 . A well-understood 
method of normalization is used to make the unsatisfactory p it j matrix easier to use. Each value along 
the main diagonal is set to 0 each line standardized. Thus, a more useful matrix, P, is built: 
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(17) 



P 



n-1,0 



P 



n-1,1 



P 



n-1,2 



0 



where i, j = 0, 1, 2, ... 

Matrix multiplication between Err t and Py brings a matrix of amended values, defined as the Err t 
matrix. The corrected matrix, C, is determined by the simple combination of the rough value matrix, 
A, and the amended value matrix, Err t , which can be simply represented by the following equation: 

C = A + Errt. 

3.4. Process of Network Error Correction 

Based on the NEC model, we can assume that the measured data of Sensor^ and Sensor j are Xi 
and Xj, respectively, which both follow Gaussian distribution. Their density functions, Pi(x) and Pj(x), 
respectively, can be used for the performance description of the sensor arrays. Xj and Xj are observed 
values of Xi and Xj, respectively. 

We introduce fiducial distance measurement, in order to reflect the size of the error between Xi and xf 



and A is the area under the density probability plot Pi(x\xj) over the interval, (x^Xj), as shown in 
Figure 3. 




j 



(18) 



where 




(19) 



Figure 3. Schematic diagram of the confidence distance measurement. 
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As mentioned above, Cor^ is defined as the fiducial distance measurement between the values of 



Sensor i and Sensor- j, where 0 < di 3 < 1. The smaller d^, the more similar the data Sensor \ and 



Sensor j are, which means that Cor^ describes the degree of fusion between Sensori and Sensor j 
Cor\j can be calculated by error function, erf (9), which is: 

erf 



Cor-ij 



(20) 



In the Learning Mode, three thresholds of Cor i3 - — namely e 2 , e and e x — can be obtained after repeated 
trials. Then: 



(21) 



1 


Po,i 


P0,2 • 


" " P0,n-1 


Pl,0 


1 


Pl,2 • 


• • Pl,n-1 


P2,0 


P2,l 


1 


• • P2,n-1 



where: 



Pn-1,0 Pn-1,1 Pn-1,2 

0, 



pij 



dij > E 2 



I ~ I (fcr ) ' £l < ^ < 5 

1, djj < £i 



(22) 



Subsequently, the standardized relationship matrix, P^, can be acquired in the 



standardization process: 



Pij 



0, 

Pij 



% = J 



(23) 



The data from one sensor are labeled as valid if the relationship value in the standardized relationship 
matrix is large enough, which means that this sensor is supported by enough other sensors. 
Otherwise, if one sensor is supported by few or no sensors, its data are labeled invalid, which means 
that the data should be deleted. The set of all the valid data is called the fusion set L = {x 1: x 2: x t }, 
where the number of elements is called the optimum fusion number, /. 

Based on the relationship matrix and fusion set, the maximum likelihood method has been adopted as 
the data fusion method. If 

1 



Pi(xi\0) 



27T<7 ; 



exp 





~Xi - 6' 


to 






J 



where % = 1, 2, I, then: 



L(x 1 ,x 2 ,...,x l ;9) = Y[p t (x t \9) 



(24) 



(25) 



t=i 



where L{x\, x 2 , xi; 9) is the maximum likelihood function based on the Maximum Likelihood 
Estimation (MLE), which means: 



d_ 
89 



L(x 1 ,x 2 ,...,xi;9) = 9 = 0 



(26) 
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and the solution is: 



9 



where 6 is the optimal fusion data of the fusion set, L = {xi, x%, ...,xi}. 
3.5. Description of the Two-Step Error Correction Process 
3.5.1. Previous Networked Error Correction Methods 



(27) 



Figure 4. Design comparison and system improvements based on the previous NEC studies. 
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(b) Present UDS method focusing on neighbor nodes. 
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As mentioned above, the focus of the study is a novel method of error correction that lends weight 
to the argument that the accuracy of the UDS system relies heavily on the fitness and robustness of the 
function involving the association between the various patterns of the error, the position of the nodes, and 
the real-time temperature. Particularly noteworthy are previous studies on the NEC which, in general, 
consist of a Learning Mode and an Operating Mode, as shown in Figure 4a. The Learning Mode trains 
the error functions of the best fit by gathering a large number of rough data for a sufficiently long 
time. Consequently, the Operating Mode will deal with the real data gathered from the sensors in an 
industrial environment, in order to reduce the detrimental effect from the inevitable interference and to 
obtain a better response within the tolerable error. Compared with the previous methods, there are some 
significant improvements in the present work, as detailed in the following Section. 

3.5.2. Design Comparison and System Improvements 

In Figure 4, it shows previous NEC systems and the present UDS method. It also illustrate the design 
comparison and the system improvements. We exploit not only the time domain model processing for 
each solo node, but also the relation matrix of neighboring nodes. Consequently, the present NEC system 
has three main advantages over the previous studies: 

• Accuracy: instead of single-ended input, the way of input the present NEC system has adopted is 
differential input, which reduces measurement error, such as cosine error and Abbe error. 

• Robustness: if one node is broken, previous NEC systems have no choice but to abandon its 
measured value, because the value is intolerant. However, the present NEC system can rationally 
estimate the value on the basis of the updated relation matrix. 

• Renewability: in previous NEC systems, once the Learning Mode is over, the function of each 
node is be irrevocable, unless starting the Learning Mode over again. However, as shown in 
Figure 4b, the relation matrix is be updated through feedback cross-talk every measurement period, 
which means that there is no boundary between the Learning Mode and the Operating Mode. 
Consequently, when the environment has changed, for example, from noon to midnight, the present 
NEC system can automatically update the relation matrix in order to adapt to the external influence. 

4. System Design of UDS 

4.1. Design of the Sensor 

The UDS is composed of three galvanically isolated main parts: a power module, a power amplifier 
module and a transmitter-receiver module. The power module provides two main functions: driving the 
sensors and ground protection. The power amplifier module, which is composed of a power amplifier 
and an inverting amplifier, can be integrated into a differential amplifier. Compared with the single-ended 
signal and the common-mode signal, the different-mode signal provides the following advantages: 

• A small signal can be more easily detected by means of controlling a reference voltage. 
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• Since one specific interference source homogeneously affects both sides of the differential signal 
to a large extent, the different-mode signal is virtually immune to electro-magnetic interference. 

• Because the different-mode system does not need to build a Virtual Ground at any point between 
ground and power when processing a bipolar signal in a single-supply system, the fidelity of the 
signal is much better than the single-ended signal and the common-mode signal. 



Figure 5. View of the design of the UDS. 
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The transmitter-receiver module is the most essential part of the application of the signal transduction. 
In the transmitter mode, the transmitters generate an ultrasonic beam; after the ultrasonic beam hits the 
target and rebounds, the addition of the two waveforms gives the return echo. In the receiver mode, the 
receiver acts as a mechanical detector of the reflected wave and outputs the echo signal already received. 
In Figure 5, it presents a view of the design of the UDS. 

4.2. Architecture of the Sensor Arrays 

Data of characteristic parameters measured by different sensors, even of the same type, are different 
to some degree. This deviation has two main causes: (1) sensor accuracy and (2) the mathematical 
algorithm adopted in data processing. The most important parameter in the UDS system is the distance 
data measured by the sensors, so we apply the same index parameter to multi- sensor measurement for 
the following reasons: 

• Using many sensors with different performance and accuracy can complement each other's 
advantages and cover disadvantages. 

• The redundancy configuration of sensors can improve reliability and measurement accuracy. 
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Figure 6. Architecture of the sensor arrays. 
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In Figure 6, it represents the architecture of the sensor arrays. Assuming that all the sensors are 
arranged as n-by-n matrices, log 2 n signals are needed for row decoding and log 2 n are for column 
decoding. Therefore, 2 log 2 n signals are all that is needed for locating any sensor at any position. There 
are two types of data the sensor arrays will transmit to the DC during data processing: the position 
signal and measured data. After cross communication between the row decoded signals and column 
decoded signals, the one-to-one relationship between the position of any sensor and its measured data 
can be established. 

4.3. Interconnection of the Sensor Arrays 

4.3.1. Schematic of the Process of Selecting An Element 

Figure 7. The equivalent reduced circuit when an element (EBA) is selected, (a) Proposed 
connection scheme with all the sensors having one end at a row line and another end at a 
column line, (b) Schematic of the proposed circuit showing the EBA filled with hashed lines. 




(a) 



(b) 



Sensors 2013, 13 



11834 



The proposed circuit for a 4 x 4 sensor array is represented in Figure 7a for illustration of the scheme. 
However, the analysis is applicable to any N x M sensor array. As shown in the Figure 7a, the columns 
and rows are connected to two digitally controlled single pole double throw switch banks, which allows 
connecting of any column to the load resistor, R L , and any row to the output node of the operational 
amplifier, with all other rows and columns staying connected to the ground. The element being accessed 
comes in the negative feedback path of the Operational Amplifier (Op Amp), as represented in Figure 7b. 
The other N — 1 sensors connected to the selected column make a parallel combination across the two 
inputs of the Op Amp, with non-inverting inputs connected to the ground. Similarly, M — 1 elements 
connected to the selected row get connected to the output node of the Op Amp with their other end 
grounded. The rest of the (N — 1) x (M — 1) elements, not physically connected to the Elements Being 
Accessed (EBA), have both of their ends at the ground [29]. 

4.3.2. Neighbor Discovery Algorithm 

The senor network can be modeled as a weighted directed graph (G), of which the nodes belong to the 
point set (V), while the network communication links belong to the edge set (E(G)). The communication 
link between Node u and Node v is defined as e(u, v), whose value of weight is defined as the data 
packet acceptance rate. According to the length of the forward path, we have defined four types of 
neighbor nodes: 

(1) The theoretical neighbor of Node v: (N(v)) : 
{u\u e V{G) A e(v, u) e E(G)}. 

(2) The non-forward neighbor of Node v: (Nq(v)) : 
{u\u G V(G) A e(v, u) G E(G) A e(u, v)} . 

(3) The one- step-forward neighbor of Node v: (Ni(v)) : 

{u\(u e N Q (v)) V(«G V(G) A e(v,u) e E(G) A 3m e V(G) A e(u,m),e(m,v) e E(G) A e(u,v))} . 

(4) The two-steps-forward neighbor of Node f:(A 2 (f)) : 

{u\(u e Ni(v)) V(«G V{G) Ae(v,u) E E(G)A 

3m, n G V(G) A e(u, m), e(m, n), e(n, v) G E{G) A e(u, v))} . 

Under the precondition that the topological structure of the wireless sensor network is known, we can 
obtain the neighborhood discovery algorithm of the UDS as shown in Algorithm 1 . 
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Algorithm 1 Procedure Calc -in- onb(G, N 0 , N u N 2 , IN 1 ,IN 2 ) 
Input: G; 

Output: N 0 ; N\, N 2 ; INi. The one-step-forwarding neighbor of Node v; IN 2 : The two-step- 
forwarding neighbor of Node v ; 
l: for all v G V{G) do 
2: N Q (y),N x (v),N 2 (y)^<l> 
3: for all e(u, v) G E(G) do 
4: if e{u,v) G £(G) then 
5: N 0 {v),N 1 {v),N 2 {v)^u 

6: else if 3m G V^G) A 3e(«, m), e(m, v) G -E(G) then 

7: iVi(v),iV 2 (V)<«-« 

8: else if 3m, n G V(G) A 3e(u,m),e(m,n),e(n,u) G -E(G) then 

9: N 2 (v) <- u 

10: end if 
11: end for 
12: end for 

13: for all v G V(G) do 
14: INi(v) <— Ni(v) — Nq{v) 
15: IN 2 {v) <- N 2 (v)-N 0 (v) 
16: end for 

5. Experiments and Results 

5.1. Comparison of Initial Data 

In order to simplify the UDS system, the experiments are based on a sample of sensor arrays, 
represented as matrices on the order of twenty by twenty. It makes the comparison of initial data between 
the ideal response, response in the absence of noise, and response with signal noise. The height of the 
semi-finished products is 50 cm, so the default distance of the parallel lines with none of the products 
passing by in the workshop that we used for experiments is also 50 cm. Therefore, the ideal value matrix 
has elements with value either 50 cm or 100 cm. Considering that interference highly depends on the 
position of the sensors in this model, we have obtained the initial data. 

Figure 8. Comparison of initial data. 




(a) Ideal response in the absence of noise. 



(b) Initial data with signal noise. 
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In Figure 8a, it shows the case of ideal response in the absence of noise, for which there are 
only two possible values, depending on which mode the sensor is in. This demonstrates that at a 
specific time, whether there is a semi-finished product passing by the sensor or not, the response 
is irregular. By comparing histograms in Figure 8a and Figure 8b, we can observe that there is 
macroscopic difference between the initial data and the ideal response, which means the measuring error 
is statistically significant. Furthermore, the mean absolute percentage error was 11.04%, which would 
not be totally acceptable in an industrial setting. Thus, in Figure 8, it shows that, in order to obtain a 
better performance, it is necessary to compensate errors resulting from the noise of the environment. 

In Figure 9, it shows the front and side views of a three-dimensional representation of the initial data 
errors, which is a better visualisation than that in Figure 8, since it is in relative terms and demonstrates 
more clearly the association between measurement error and noise from the environment. Furthermore, 
in Figure 9, it also shows that there is little difference between the front view and the side view. One 
plausible explanation is that the error function graph of each node is approximately symmetrical with 
rotation around the perpendicular. Based on these factors, we decided to focus on the front view for 
this model. 

Figure 9. Initial data errors. 
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(a) Three-dimensional diagram. (b) Front view. (c) Side view. 

5.2. Analysis of the Relationship between Nodes in Not Z-MODE and Nodes in Z-MODE 

From Figure 10, we can observe that the correlation coefficient fluctuates a great deal with the 
change of the relative position between the node in Not Z-MODE and the nodes in Z-MODE. Although 
this observation cannot establish cause and effect, it strongly suggests that the size of interference is 
largely related to relative position. Additionally, one of the plausible interpretations of this observation, 
considering the field of a manufacturing workshop shown in Figure 1, is that the whole production 
process can be decomposed into several smaller, repeated processes, where interference from the 
environment in subprocess / is approximately parallel to subprocess J. This observation lends significant 
weight to the argument that it is reasonable for the value of an arbitrary node to be modified through a 
two-step process: 

(1) Select a set of data for which the correlation coefficients are large enough. 

(2) Correct this value of the arbitrary node based on the set data chosen in the first step. 
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Figure 10. Correlation coefficient between a node in Not Z-MODE and the nodes in Z- 
MODE, obtained in the Learning Mode. 




5.3. Process of Error Correction 

In Figure 11, it shows that the experimental response is more accurate than the initial data 
shown in Figure 9. Although some nodes' values were close to zero, the others were undesirably 
large. Furthermore, the reason for this observation is that we have only modified the nodes in Not 
Z-MODE and not the nodes in Z-MODE. Consequently, additional optimization is required to obtain a 
desirable response. 

Figure 11. First-time corrected error values by applying a UDS. 
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In Figure 12, it shows that the experimental data have been fitted to ideally respond after modifying the 
nodes in Not Z-MODE. Compared with Figure 9b, the optimized method of distance measurement with 
the capability of NEC performs excellently, since it takes full advantage of the correlation coefficient 
among neighbor sensors. So far, we have described the process of one application to the UDS. However, 
determining a reasonable number of iterations is of great importance. 
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Figure 12. Second-time corrected values by applying a UDS. 
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5.4. Optimal Number of Iterations 

In Figure 13, it shows the tendency of the mean absolute percentage error along with the increasing of 
iterations, from which we can observe that the mean absolute percentage error has significantly reduced, 
from 0.08% to 0.01% in the first three iterations. This observation demonstrates that the adaptive NEC 
system has the ability to eliminate the steady state error. Moreover, the mean absolute percentage error 
uniformly converges to zero once the number of iterations is large enough, and the convergence rate is 
approximately O(ln). The reduced rate suddenly slows down from the fourth iteration, however. Taken 
together, this indicates that the most reasonable number of iterations is three, since this gives a desirable 
response without extra cost. 

Figure 13. Mean absolute percentage error along with the increasing of iterations. 
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6. Conclusions 

A novel method based on NEC for obtaining a satisfactory distance measurement has been proposed 
in this work. It can provide a high-speed and high-accuracy solution for distance measurement, which 
is vital in most industrial fields. One of the most obvious flaws of the previous methods is that none of 
them has considered the indirect association among neighboring sensors, which leads to less accuracy 
and extra cost. Therefore, one of the main aims of this study was to propose a novel method based on 
NEC that can provide a high-speed and high-accuracy solution for distance measurement. We applied 
the maximum likelihood method to determine the optimal fusion data set and used a neighbor discovery 
algorithm to identify neighboring nodes at high speed; furthermore, we adopted the NEC optimization 
algorithm, which takes full advantage of the correlation coefficient among neighboring sensors. The 
experimental results demonstrate that the ranging errors of the NEC system are within 2.20%, and the 
mean absolute percentage error is reduced to 0.01% after three iterations of this method. 

Our proposed method has several strengths. However, several limitations merit comment. The most 
critical of them is that we are unable to establish and then test the exact cause and effect equation of 
the numerical error and the propagation of the ultrasonic wave affected by the environment, such as 
temperature, humidity and atmospheric pressure. Nevertheless, in spite of this limitation, we believe 
this study opens new paths of investigation in distance measurement using a UDS and, to some 
extent, modifies the manner in which we understand the method of error compensation, because new 
insights into it will likely generate novel methods to measure distance more effectively and quickly than 
previously. There is, therefore, a great need for further research in this area. 

Acknowledgments 

This work has been partially supported by the project "Science and Technology Plan of Shandong 
Province, China (No.2012GB020108)" and the project "National Science Foundation of China 
(No. 61070022)". 

Conflicts of Interest 

The authors declare no conflicts of interest. We declare that we have no financial and personal 
relationships with other people or organizations that can inappropriately influence our work, and there is 
no professional or other personal interest of any nature or kind in any product, service and/or company 
that could be construed as influencing the position presented in or the review of the manuscript. 

References 

1. Escol, A.; Planas, S.; Rosell, J.R.; Pomar, J.; Camp, F.; Solanelles, F.; Gracia, F; Llorens, J.; Gil, E. 
Performance of an ultrasonic ranging sensor in apple tree canopies. Sensors 2011, 11, 2459-2477. 

2. Andjar, D.; Weis, M.; Gerhards, R. An ultrasonic system for weed detection in cereal crops. 
Sensors 2012, 12, 17343-17357. 

3. Majchrzak, J.; Michalski, M.; Wiczynski, G. Distance estimation with a long-range ultrasonic 
sensor system. IEEE Sens. J. 2009, 9, 767-773. 



Sensors 2013, 13 



11840 



4. Giannoccaro, N.; Spedicato, L.; di Castri, C. A new strategy for spatial reconstruction of orthogonal 
planes using a rotating array of ultrasonic sensors. IEEE Sens. J. 2012, 12, 1307-1316. 

5. Lee, J.R.; Chia, C.C.; Shin, H.J.; Park, C.Y.; Yoon, D.J. Laser ultrasonic propagation imaging 
method in the frequency domain based on wavelet transformation. Opt. Lasers Eng. 2011, 49, 
167-175. 

6. Wang, X.; Tang, Z. Optional optimization algorithms for time-of-flight system. IEEE Trans. 
Instrum. Meas. 2011, 60, 3326-3333. 

7. Zhao, H.; Song, P.; Urban, M.W.; Kinnick, R.R.; Yin, M.; Greenleaf, J.F.; Chen, S. Bias observed in 
time-of-flight shear wave speed measurements using radiation force of a focused ultrasound beam. 
Ultrasound Med. Biol. 2011,37, 1884-1892. 

8. Canali, C; de Cicco, G.; Morten, B.; Prudenziati, M.; Taroni, A. A temperature compensated 
ultrasonic sensor operating in air for distance and proximity measurements. IEEE Trans. Ind. 
Electron. 1982, IE-29, 336-341. 

9. Carullo, A.; Ferraris, E; Graziani, S.; Grimaldi, U.; Parvis, M. Ultrasonic distance sensor 
improvement using a two-level neural-network. IEEE Trans. Instrum. Meas. 1996, 45, 677-682. 

10. Angrisani, L.; Baccigalupi, A.; Moriello, R.S.L. A measurement method based on Kalman filtering 
for ultrasonic time-of-flight estimation. IEEE Trans. Instrum. Meas. 2006, 55, 442^448. 

1 1 . Jiang, S.B.; Yang, C.M.; Huang, R.S.; Fang, C.Y; Yeh, T.L. An innovative ultrasonic time-of-flight 
measurement method using peak time sequences of different frequencies: Part I. IEEE Trans. 
Instrum. Meas. 2011, 60, 735-744. 

12. Hajiyev, C. Innovation approach based measurement error self-correction in dynamic systems. 
Measurement 2006, 39, 585-593. 

13. Yan, T.H.; Wang, W.; Chen, X.D.; Li, Q.; Xu, C. Design of a smart ultrasonic transducer for 
interconnecting machine applications. Sensors 2009, 9, 4986-5000. 

14. Phan, PA.; Gale, T. Two-mode adaptive fuzzy control with approximation error estimator. IEEE 
Trans. Fuzzy Syst. 2007, 15, 943-955. 

15. Gu, L.; Kou, X.; Jia, J. Distance Measurement for Tower Crane Obstacle Based on Multi-Ultrasonic 
Sensors. In Proceedings of 2012 IEEE International Conference on Automation Science and 
Engineering (CASE), Seoul, Korea, 20-24 August 2012; pp. 1028-1032. 

16. Tsai, T.M.; Yen, PH. Improvement in stage measuring technique of the ultrasonic sensor gauge. 
Measurement 2012, 45, 1735-1741. 

17. Srinivasan, S.; Pandharipande, A. Self-configuring scheduling protocol for ultrasonic sensor 
systems. IEEE Sens. J. 2013, 13, 2517-2518. 

18. Kumar, R. Game theoretic pattern analysis for identification of odors/gases using response of a 
poorly selective sensor array. IEEE Sens. J. 2013, 13, 1110-1116. 

19. Lee, K.Y; Huang, C.E; Huang, S.S.; Huang, K.N.; Young, M.S. A high-resolution ultrasonic 
distance measurement system using vernier caliper phase meter. IEEE Trans. Instrum. Meas. 
2012,67,2924-2931. 

20. Bischoff, O.; Heidmann, N.; Rust, J.; Paul, S. Design and implementation of an ultrasonic 
localization system for wireless sensor networks using angle-of-arrival and distance measurement. 
Proced. Eng. 2012, 47, 953-956. 



Sensors 2013, 13 



11841 



21. Zhang, X.B.; Yang, B.C.; Zhang, T.; Zhang, S.Y. Local influence on the error-correction variable 
in a cointegrated system. J. Syst. Eng. Electron. 2001, 12, 1-8. 

22. Singh, M.; Kumar, P. An efficient forward error correction scheme for wireless sensor network. 
Proced. Technol. 2012, 4, 737-742. 

23. Luo, J.; Shi, C; Shi, W.; Ling, Q. Joint Node Localization and Error Ccorrection in Wireless 
Sensor Networks. In Proceedings of the 31st Chinese Control Conference (CCC), Anhui, China, 
25-27 July 2012; pp. 6600-6604. 

24. Yu, K.; Barac, R; Gidlund, M.; Akerberg, J. Adaptive Forward Error Correction for Best 
Effort Wireless Sensor Networks. In Proceedings of 2012 IEEE International Conference on 
Communications (ICC), Ottawa, Canada, 10-15 June 2012; pp. 7104-7109. 

25. Barcelo-Lladd, J.E.; Perez, A.M.; Seco-Granados, G. Enhanced correlation estimators for 
distributed source coding in large wireless sensor networks. IEEE Sens. J. 2012, 12, 2799-2806. 

26. Zhang, J.; Shen, X.; Zeng, H.; Dai, G.; Bo, C; Chen, E; Lv, C. Energy-efficient and localized lossy 
data aggregation in asynchronous sensor networks. Int. J. Commun. Syst. 2012, 26, 989-1010. 

27. Malek, N.I.; Ijardar, S.P; Master, Z.R.; Oswal, S.B. Temperature dependence of densities, speeds 
of sound, and derived properties of cyclohexylamine + cyclohexane or benzene in the temperature 
range (293.15 to 323.15 K). Thermochim. Acta 2012, 547, 106-119. 

28. Tang, J.; Liu, J. Variance of speed of sound and correlation with acoustic impedance in canine 
corneas. Ultrasound Med. Biol. 2011,57, 1714-1721. 

29. Saxena, R.; Bhan, R.; Aggrawal, A. A new discrete circuit for readout of resistive sensor arrays. 
Sens. Actuators A: Phys. 2009, 149, 93-99. 

© 2013 by the authors; licensee MDPI, Basel, Switzerland. This article is an open access article 
distributed under the terms and conditions of the Creative Commons Attribution license 
(http://creativecommons.Org/licenses/by/3.0/). 



